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1. Introduction 

This work is motivated by the paper [LPj . In this paper, the au- 
thors consider the scaling limits of three discrete aggregation models 
with multiple sources - the internal DLA (diffusion-limited aggrega- 
tion), the rotor- router, and the divisible sandpile model. They show 
that for all three models, the scaling limit is the same, and it is a so- 
lution of the Hele-Shaw injection problem with multiple sources. In 
two dimensions, the latter problem can be solved explicitly using the 
theory of quadrature domains, which gives an explicit formula for the 
scaling limit. 

They also consider the same aggregation models in two dimensions 
with just one source, but with an additional condition on the posi- 
tive (horizontal) half-axis (namely, the condition that a particle hit- 
ting the positive half-axis is killed, or the condition that it is directed 
downward). In these cases, the existence of the scaling limit remains 
unknown, although it is expected to exist, and computer-generated 
pictures for its shape are given in Fig. 4 of [LPj . 

The goal of this paper is to study the Hele-Shaw problems which are 
expected to be the scaling limits of the above models with a condition. 
Namely, we provide an explicit solution for the Hele-Shaw problem 
corresponding to the killing condition, which is a close fit with the 
left shape at Fig. 4 of [LPj . We also describe moment properties of 
the solution of the Hele-Shaw problem corresponding to the downward 
condition (the right shape at Fig. 4 of [LPj ). although we are unable 
to compute this shape explicitly. 
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2. Flows with the killing/reflecting condition 

2.1. The killing condition on a half-axis. Consider a 2-dimensional 
discrete aggregation model with injection at the origin and killing con- 
dition on a half-axis; we prefer to impose this condition on the negative 
half-axis (rather than the positive half-axis used in [LP] ). To be more 
concrete, consider the internal DLA model. In this model, particles 
are emitted from the origin into a rectangular lattice, and each parti- 
cle walks randomly on this lattice until it reaches an unoccupied node, 
where is stays, and a new particle is emitted. In addition, one stipulates 
that a particle that hits the negative half-axis is killed (after which a 
new particle is emitted). The scaling limit of this model is the limit in 
which the number N of surviving particles Q goes to infinity, and the 
spacing of the lattice goes to zero as (A/N) 1 / 2 , where A is a constant. 
If the scaling limit exists, it is a region Q in the plane of area A, con- 
taining the origin in its interior. We would like to compute this region 
explicitly using conformal mappings, assuming the scaling limit exists 
and is the solution of the Hele-Shaw problem naturally associated to 
this model. 

Now consider the Hele-Shaw problem for the asymptotic region Q. 
We refer the reader to the books |VEj and [GVJ for the basics on Hele- 
Shaw moving boundary problems. 

Problem 2.1. Let $ be the harmonic potential in Q with deleted 
negative half-axis, with boundary conditions $ = on the boundary 
r of Q and also on the negative half-axis (which is the continuous 
counterpart of the killing condition on the negative half-axis in the 
discrete model), and $ ~ — Re(,z -1 / 2 ) near zero (where z is the complex 
coordinate in the plane). The condition on Q is that the infinitesimal 
deformation of the boundary given by the normal derivative of $ is 
homothetic. 

Let us find fl from this data (up to scaling). For this purpose, let 
/ : D — > Q be the conformal map of the unit disk onto Q with /(0) = 
and f'(Q) > (it is unique). Let = $(/(«)). Then * is a 

1 Note that most of the particles will not survive. Namely, as was pointed out to 
us by Lionel Levine, it follows from the Bcurling estimate that for N particles to 
survive, the number of produced particles will have to be ~ CN 5 / 4 , where C is a 
constant. 



ON HELE-SHAW PROBLEMS ARISING AS SCALING LIMITS 



3 



harmonic potential on the disk with a cut along the negative half-axis, 
and \I/ behaves like — Re(w -1 / 2 ) at zero. Thus, up to scaling, 

(we choose the standard branch of the square root with a cut on the 
negative half-axis). Indeed, \1/ is uniquely determined by its properties, 
and the given \1/ satisfies them. 

Now, if |C| = 1 then the velocity vector of the boundary at the 
point f(() under the transformation by the normal derivative of $ is 
ff l/'(C)l _1 ) while the outward unit normal to the boundary at f(() is 
tjttqi ■ Thus the condition of homothetic deformation has the form 



i.e., 



(1) Re(c/'(C)/(C)) = ^,ICI = i, 

or 

(2) Re(C/'(C)7(0) = \(C 1/2 + C 1/2 ), ICI = l. 

Now we would like to solve (j2J). We use the method similar to one 
described in [EEK] . for evolution of polygons (Section 9). 

Let D = (d{. Also for any function g set g*(() = g(l/(), so (Dg)* = 
-D{g*). 

Note that / is real, i.e. it has real Taylor coefficients at 0. Thus 
equation (j2J) can be rewritten in the form 

(3) Df-r+f. (Dfy = c i/2 + c 1/2 , ici = i. 

Differentiating this, we get 

(4 ) D 2 f-r-f-(D*fy = l(c /2 -c l/2 ). 

Let 

h = 7 ' 

Then equations fl3El) imply that 

h-h* = o,\C\ = i. 

But by its definition, h extends to a function which is holomorphic 
everywhere in the unit disk except possibly the point —1. Hence h is 
holomorphic in the Riemann sphere with the possible exception of the 
point —1. Also, it's clear that Df(—1) = 0, which by the removable 



4 



PAVEL ETINGOF 



singularity theorem and Liouville's theorem implies that in fact h is a 
constant. 

Thus, we find that / satisfies the differential equation 

D 2 f + - ] —^Df = hf 
J 21 + C 



The constant h is easy to find from the condition that f(() is propor- 
tional to ( + 0(( 2 ): namely, h = |. 

From this is it easy to solve the differential equation by the power 
series method. The answer is, up to scaling 

fx) = Uly (-0*. 

AW 16^ (k 2 - 1/4) (A; 2 - 9/4) 1 ^ 

(the normalization is such that /'(0) = 1). Thus, we have 

Proposition 2.2. The region solving Problem \2.1\ is defined by the 
conformal mapping 

(5) /(C) = |(1 + CHC 1 - (1 - OC^arctanC 1 / 2 ) - ~. 

Thus we have a "logarithmic cusp" near ( = —1, i.e. /(— 1 + u) 
behaves like cu 2 \ogu (as opposed to the ordinary semicubic cusp, with 
local behavior cu 2 ). We also see that the region is twice as thick in the 
positive direction as it is in the negative direction. 

Formula (151) is a very close fit with the left shape on Fig. 4 of |LPj 
(when it is rotated by 180° to match the positive and the negative half- 
axes). This leads to a conjecture that the region defined by formula 
(j3j) is the scaling limit of the model. 

2.2. The killing condition on the sides of an angle. The analysis 
of the previous subsection can be generalized to the following more gen- 
eral problem. Let < b < 1, and consider the angle {z e C|argz e [— 7r6, i 
Consider the internal DLA model as above inside the angle with a 
killing condition on the sides of the angle. Then the setting of the 
previous subsection is the special case 6=1. 

The Hele-Shaw problem for the asymptotic region Q in this case is 
as follows. 

Problem 2.3. Suppose that Q is a wedge-shaped region bounded by 
the sides of the angle and some curve T connecting them (symmetric 
with respect to the horizontal axis). Let $ be the harmonic potential 
in Q which vanishes on the boundary T and the sides of the angle, and 
behaves like — ~Re(z~ l l 2h ) near zero. Then the defining condition for Vt 
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is that T transforms homothetically under the infinitesimal deformation 
by the normal derivative of $. 

Remark 2.4. Note that the condition 6 < 1 is not essential, as for 
6 > 1 we can put the angle (which is > 2n) on the helical Riemann 
surface, which is the universal covering of the complex plane punctured 
at zero. 

The region Q for any b can be found similarly to the case 6=1 
considered above. Namely, let / be the conformal map of the circular 
sector defined by the inequalities \z\ < 1, — nb < argz < nb (lying on 
the Riemann surface if b > 1) onto Q, such that / bijectively maps 
straight sides to straight sides (so /(0) = 0), and f'(0) > 0. Then we 
have an expansion 

/(0 = 5>c 1+ * /6 , 

k>0 

for some at G M. 

By a method analogous to the case b = 1, we derive a differential 
equation for /, which has the form 

1 1 - C l l b 

D 2 f + —rjrDf = (1 + 1/26)/. 

26 1 + C l/ 

Then we can solve for / by the power series method, and up to scaling 

we get 

Proposition 2.5. The region solving Problem \2.3\ is defined by the 
conformal mapping 

(6) f(0 = (F( aj [3,r,-( 1/b ), 

where F is the Gauss hypergeometric function: 

n>0 j=0 K J J \ i J / 

and 

a = -1/2, p = 2b, 7 = 26 + 3/2. 

Example 2.6. Consider the special case 6 = 1/2 (the angle is n, i.e. 
its boundary is the imaginary axis). In this case the solution is 

/(C) = Cn-l/2,l,5/2;-C 2 ), 

i.e., 

/(C) = |((C - C 1 ) + (C + CTarctanC). 
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Remark 2.7. This analysis is similar to the analysis in [GVJ, subsec- 
tion 2.2.3, where solutions are also expressed via the Gauss hypergeo- 
metric function. 

Lionel Levine has simulated on a computer the corresponding dis- 
crete (rotor-router) model for 6=1/2 and 6 = 1/4 (the dynamics in the 
half-plane and the quarter-plane, respectively), and got a very close fit 
with formula (J6j). This gives rise to a conjecture that for these values 
of b (and perhaps for all) the scaling limit is given by formula (jH]). 

2.3. The killing and reflecting case. Another problem one can 
consider is the same as the previous one, except that the angle is 
{z G C|argz G [0, Tib}} (i.e., the upper half of the angle considered be- 
fore), and on the right side of the angle (the positive half-axis) we have 
not the Dirichlet boundary condition for $, but the Neumann boundary 
condition, = 0. This can be conjecturedly interpreted as the scaling 
limit of the random walk problem on a rectangular lattice in this angle 
with a killing condition on the left side and reflecting condition on the 
right side. 

The problem of finding Q in this case reduces to the previous one, 
by doubling the angle using reflection with respect to the positive half- 
axis. In particular, if b = 2 (the angle is the full plane, and we have 
the reflecting condition above the positive half-axis and the killing con- 
dition below), then the conformal map of the disk with a cut at the 
positive half-axis to the region Q is 

/(C) = CF(-l/2,4, 11/2; -C 1/2 ), 

where C 1/2 := r l / 2 e ie ' 2 , 6 G [0,2tt). 

3. Flows with the killing-passing-reflecting condition 

Consider now the discrete model as above, with the condition on the 
positive half-axis saying that particles coming from above are reflected, 
while particles coming from below pass through with probability p and 
are killed with probability 1 — p. If p = 0, this is the killing-reflecting 
case considered in the previous section, and if p = 1, this is the passing- 
reflecting case, whose simulation is shown on the right side of Fig. 4 of 

We expect that the scaling limit of this process exists for any p G 
[0, 1], and is described as follows. 

Problem 3.1. Let Q be the asymptotic region. Then the boundary 
of Q consists of a curve T and the interval [a, b] of the positive half- 
axis (the curve T goes counterclockwise from 6 to a around the origin). 
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Then the boundary conditions for the potential <3> in Q are: 



on r; 



($„)+ = on [a, &]; 
$_ = on [0,a]; 
($„)+=p($ y )_ on [0,a]. 



Here is the derivative of $ with respect to y, and +, — denote the 
one-sided limits from above and below (<& and $ y are not assumed 
continuous on (0, a)). 

The condition on Q, as before, is the homothetic transformation 
under the flow defined by the normal derivative of $. 

Remark 3.2. 1. To motivate these boundary conditions, it is conve- 
nient to use the following reformulation of the killing-passing-reflecting 
condition in the discrete model: all the particles hitting the positive 
half-axis from below are killed, but for each k particles killed on a small 
interval / we produce about pk particles which originate at / and move 
in the upward direction. This explains the second and third boundary 
conditions (reflecting condition on [a, b] from above, killing condition 
on [0, a] from below), and clarifies the meaning of the fourth boundary 
condition, which says that the flux through [0, a] on the upper side is 
p times the flux on the lower side. 

2. Note that for p — 0, Problem 13.11 reduces to the problem of 
Subsection 12.31 for 6 = 2. 

Unfortunately, we did not manage to compute the solution Q of Prob- 
lem 13.11 explicitly. However, here is an interesting moment property of 
this solution (which we expect to determine fl uniquely up to scaling). 

Define a continuous branch of log(z) in Q \ M + by defining it to be 
real on the reflecting boundary (i.e. the limit of log(z) from above is 
real for z G 7^ 0). This also defines branches of z s = e slog ^ for 

any s. 

Proposition 3.3. (%) If Q is a solution of Problem \3.1l then one has 




where a = arccos(p)/27r. In particular, < a < 1/4. 
(ii) For p = 1 (a = 0), one has 
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Proof, (i) Let u be a function on Q \ M.+. Then the time derivative of 
the moment f n udxdy under the Hele-Shaw flow is given by 

— I udxdy = f u——dl = 
dt J n J r dn 

r <9$ r b r a 

/ u——dl — / u + (§ y ) + dx + / U-(§ y )_dx. 
Jan dn J ' J 

By Green's formula, this equals 



[u 

n 



f du f b [ a 
A§-§Au)dxdy + / $— dl - / u + (® y )+dx + / U-{§ y )-dx. 

Jan v* 1 Jo ' Jo 



If u is harmonic and vanishes to sufficient order at then the first 
summand is zero. So, using the first three boundary conditions, we 
have 

d f f b f a f a 

(7) — / udxdy = / <& + {u y ) + dx— I u + (§ y ) + dx+ / 
dt Jq Jo Jo Jo 

Now suppose that u = Ke(z n±a ) = r n+a cos((n ± a)9), z = re 10 , 9 G 
[0, 27r). Then (u y )+ = 0, so the first summand in ([7j) vanishes. Thus, 
we have 

d f f a f a 

(8) — / udxdy = / U-($ y )-dx — / u + (& y ) + dx + . 
dt Jn Jo Jo 

But we have u + = x a , m__ = x a cos(27ra) = px a . So by the fourth 
boundary condition, we get from (jTJ) that ^ f Q udxdy = 0. But since 
the time evolution of Q is homothetic, and J n udxdy is homogeneous 
of nonzero degree under dilations, we get J n udxdy = 0, as desired. 

(ii) follows from (i) by considering the difference of the two equations 
in (i) and computing the first term of the Taylor expansion in a at 
a = 0. □ 

Remark 3.4. Let Q be a region solving Problem 13. II for p = 1 (which 
is expected to match the right shape at Fig. 4 of [LP]), and Q be the 
complex conjugate region. Proposition I3.3( ii) implies that 

z n dxdy + / z n dxdy = 0, n E N. 
n Jn 

This implies that the "balayage sum" E of Q and Q (i.e., the region 
whose moments are equal to the sum of the moments of Q and of Q, 
see |GVj ) satisfies the condition J E z n dxdy = 0, n G N, i.e. is a disk 
(see |VE] and references therein). 
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Note that for p = (a = 1/4), i.e. the killing-reflecting case dis- 
cussed in Subsection [531 the two arithmetic progressions in Proposition 
I3.3( i) combine into a single one, and yield 



Jn 

This can be generalized to the case of any angle (i.e., any b). Namely, 
we have the following proposition. 

Proposition 3.5. Let Q be a region solving Problem \2.'d Then 



Jn 

Therefore, the solution of the killing-reflecting problem with parameter 
b (which is the upper half of the region Q) satisfies the equality 



The proof of this proposition is similar to the proof of Proposition 
13.31 using Green's formula. 

A similar property was established by Lionel Levine for the discrete 
model (for 6=1). 
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